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We investigate the correlated predictions of singlet extended MSSM models for direct detection of 



the lightest neutralino with its cosmological relic density. To illustrate the general effects of the 
• Jh . 

singlet, we take heavy sleptons and squarks. We apply LEP, (g — 2)^ and perturbativity constraints. 

H \ 

We find that the WMAP upper bound on the cold dark matter density limits much of the parameter 
space to regions where the lightest neutralino can be discovered in recoil experiments. The results 
for the NMSSM and UMSSM are typically similar to the MSSM since their light neutralinos have 
similar compositions and masses. In the n/sMSSM the neutralino is often very light and its recoil 
detection is within the reach of the CDMS II experiment. In general, most points in the parameter 
spaces of the singlet models we consider are accessible to the WARP experiment. 
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I. INTRODUCTION 



The evidence for cold, non-baryonic dark matter is one of the strongest cases for physics 
beyond the Standard Model (SM). The most plausible candidates are axions and weakly 
interacting massive particles (WIMPS). Many models have been proposed that provide a 

WIMP candidate for dark matter (DM}.__These models include Supersymmetry ii b, a, 



Extra Dimensions 



IB, 



Little Higgs 



and twin Higgs models [9 , llO|] . It is anticipated 



that the LHC and ILC may directly produce the DM particle, providing crucial information 
on its existence, origin and properties [ll|, . Recoils of nuclei from 



the scattering of WIMPS can also provide direct evidence for this DM candidate 



20]. 



Alternative avenues exist to discover the DM particle indirectly. These inc 
of gamma ray lines in the galactic halo peaked at the neutralino mass 



of high energy neutrinos from WIMP annihilation in the sun 23 . 



24 



detection 



26 



21, 



28 



21 

3 



ude observation 
221 ] . observation 
and antiparticle 



29]. 



The Minimal Supersymmetric Standard Model (MSSM) with conserved -R-parity has a 
viable dark matter candidate, the lightest neutralino, that can naturally explain the abun- 
dance of cold dark matter (CDM) in the universe. Other motivations for the MSSM include 
solutions to the gauge hierarchy problem, the quadratic divergence in the Higgs boson mass, 
and gauge coupling unification 30, |31_|, |32| . The MSSM lagrangian contains a Higgsino mix- 
ing parameter, /z, the only massive parameter of the model that conserves supersymmetry. 
The /i parameter must be at the electroweak (EW) or TeV scale to explain electroweak 
symmetry breaking. This is the so-called hierarchy problem of supersymmetry because a 
priori the value of /i is arbitrary 33J . The problem may be resolved by promoting the [i pa- 
rameter to a dynamical field whose vacuum expectation value (S) and coupling A determine 
the effective /i-parameter, 

/Xeff = MS). (1) 



Singlet extended mode 



sectors [34], 



35 



36 



37|, 



sfxMSSM) 



38 



39 
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have significant consequences for the Higgs and neutralino 



41 



42 



43 
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45]. 



In this paper we focus on the some of the proposed singlet extensions of the MSSM: (i) 



et term 



the Next-to-Minimal Supersymmetric SM (NMSSM) 46|, |47J, |48J with a trilinear sing 
in the superpotential, (ii) the Nearly-Minimal Supersymmetric SM (nMSSM) 49j, 150 
with a tadpole term in the superpotential, and (iii) the U{1)' '-extended MSSM (UMSSM) 



51 



52] 
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53|, |54j, |55j with an extra Z' gauge boson, as detailed in Table[T]along with the respective sym- 
metries and numbers of neutralino, chargino and Higgs states. A Secluded [/(l)'-extended 
MSSM (sMSSM) [35|, |56| with three singlets in addition to the standard UMSSM Higgs 
singlet is equivalent to the nMSSM in the limit that the additional singlet vevs are large, 
and the trilinear singlet coupling, X s , is small [37]. Thus, we collectively refer to the nMSSM 
and sMSSM as the n/sMSSM. 

Other singlet extensions are possible, and often are derived from string constructions. The 
UMSSM is an example of such string constructions 



and articles, see, e.g., Refs 



57j, 



58, 



59, 



60 



61 



62 



63 



54, 



57 



58 



59 



60] . For recent reviews 



64| . For a review of supersymmetric 



singlet models, see Ref. [34j. The SM may also be extended to include a singlet. Depending 



on the parity of the singlet, it may mix with the Higgs boson 65|, [66] or provide a viable 



dark matter candidate 



67 



68|. 



Model: 


MSSM NMSSM 


nMSSM 


UMSSM 


sMSSM 


Symmetry: 




Z 3 




u(iy 


u(iy 


Extra 




§S 3 


t F S 




\sS1S2S3 


superpotential term 




(cubic) 


(tadpole) 




(trilinear secluded) 


X° 


4 


5 


5 


6 


9 
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3 


3 
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4 


1 


2 
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1 


4 



TABLE I: Symmetries associated with each model and their respective terms in the superpotential; 
the number of states in the neutralino and Higgs sectors are also given. All models have two 
charginos, Xi > an d one charged Higgs boson, H^. 



The extended MSSM models (xMSSM) have extra neutra! 
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ino and Higgs states, that 
, and significantly influence 



affect the spectra and composition of the particle states 
the parameter space for the observation via direct or indirect detection of neutralino dark 
matter. For example, the processes involved in indirect detection have been shown to be 
radiatively enhanced by the extended Higgs sector in the NMSSM 69(. The direct detection 



prospects have been recently investigated in the NMSSM 70j, |7l). The relic density of the 



lightest neutralino can also be dependent on the model; the density relative to the closure 



3 



density is very roughly given by the total annihilation cross-section by 72] 

W - ^4 (2) 

where v is the relative velocity and a ann includes effects from neutralino co-annihilation with 
particles of similar mass. The non-baryonic Dark Matter relic density relative to the critical 
density, VL DM) is determined by the WMAP 3-year CMB data and the spatial distribution 
of galaxies to be [fjl, Q 

tt DM h 2 = 0.111 ±0.006 (3) 

where h = 0.74 ± 0.03 is the Hubble constant. However, there could be multiple origins of 
dark matter, such as neutralinos and axions or non-thermal production, so strictly speaking 
Eq. ([3]) provides only an upper bound on neutralino dark matter. 

In this paper we consider the neutralino and Higgs states of the singlet extended super- 
symmetric models. The predictions for the neutralino relic density and neutralino-nucleon 
scattering cross-section are often correlated. We find that the singlet and singlino admix- 
tures of the Higgs bosons and neutralinos, respectively, can dramatically alter the predicted 
relic density and the direct detection scattering cross-section. We calculate the xMSSM 
predictions, along with those of the MSSM 1 , with the parameterization of the Higgs and 
neutralino sectors given in Ref. 3^]. Many studies on the relic density and direct detection 
of the lightest neutralino have been performed in the CMSSM which imposes scalar mass, 
gaugino mass and soft trilinear mass unification at the GUT scale, thereby reducing the 



number of parameters in the model [7j. 



9, h y, 



771 ] . Additional work on the connection 



between the Higgs sector and astrophysical constraints in the CMSSM have also been done 
in ttie^MSSM Q3- 

More general recent analyses of the MSSM can be found in Refs. 



79|. 



In Section [TTJ, we analyze the neutralino annihilation and coannihilations and how the 
current value of flDAih 2 constrains the Higgs and neutralino spectra. In Section HTT1 we 
consider neutralino-nucleon scattering and the prospects for direct detection of neutralino 
dark matter in the xMSSM models. Finally, in Section IIVt we summarize our conclusions 
and provide an outlook on discovery for the models we consider. We reserve the Appendix 



1 Note that the version of the MSSM we adopt is more general than mSUGRA since we set parameters at 
the TeV scale. 
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for discussions of the neutralino mass matrix of the extended models and couplings which 
are altered or completely new relative to the MSSM. Further, we discuss the constraints 
from the (g — 2) M measurement and perturbativity of the couplings. 

II. NEUTRALINO ANNIHILATION AND COANNIHILATION 

The light neutralino states in singlet extended models often have an appreciable singlino 
component or are dominantly singlino. The singlino admixture has a strong influence on the 



neutralino annihilation cross-section and consequently on the neutralino relic density [381.180] . 
Coannihilations with supersymmetric particles of similar mass may also significantly affect 
the relic density. 

We calculate the relic abundance of the lightest neutralino in each model using appropri- 



ate modifications of the DarkSUSY (DS) code 8JJ 2 . We change the Higgs and neutralino 
couplings in DS to account for the additional interactions of the singlets and singlinos (and 
Z' inos in the UMSSM). We list the couplings in AppendixEl The number of Higgs and neu- 
tralino states in DS are increased to account for the additional states in the extended models 
relative to the MSSM. The Higgs boson width is recalculated for the different spectra and 
changes in couplings 3 . Our calculation of the thermally averaged annihilation cross-section 
at tree level in our modified version of DS includes all possible initial, final, and exchanged 



particles of the xMSSM 4 . For general discussion of the relic density calculation, see Ref. [84]. 

Due to coannihilations with the extra neutralino states in xMSSM models, the relic den- 
sity can be different from the MSSM result. While Boltzmann suppression makes contribu- 
tions from particles with masses > 1.4 m x o small, we include coannihilations from particles 
with masses up to 2.5 m x o (Ref. [85] included coannihilations up to 2.1 m x o). 

We fix the slepton and squark masses at a high value (5 TeV) to focus our attention 
on the effects the singlet has on the relic density. Therefore, these scalar sparticles do not 
coannihilate with the neutralino, making their effect on the relic density negligible 5 . The 



2 Other supersymmetric dark matter codes exist, including MicrOMEGAS [82] and ISATOOLS 83|. 

3 We also include the Higgs decay to off shell gauge bosons, which can be important for intermediate Higgs 
masses (120 GeV < M H% < 2M Z ). 

Note that we do not include exotic states that can be present in some extended models. For a list of the 



4 



initial, final and exchanged states considered by DarkSUSY in the MSSM, see Table 3 in Ref. 81 1. 
5 Of course, additional coannihilation solutions would be possible in all models for lighter squarks and 
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FIG. 1: Neutralino relic density versus the lightest neutralino mass. The relic density is constrained 
to be in the region 0.123 > QoAih 2 > 0.099 provided that the model is solely responsible for the 
observed dark matter. The efficient annihilations through the Higgs boson pole in the MSSM, 
NMSSM and UMSSM are evident at m x o ~ M H J2 ~ 60 GeV and of the Z boson pole at 
m x o ~ M z /2 in the n/sMSSM. 



n/sMSSM supersymmetric spectrum is typically heavy compared to the lightest neutralino 
making the contribution to the relic density from coannihilations negligible. In the MSSM, 
NMSSM, and UMSSM, if the relic density without the effect of coannihilations is at or 
above the observed value of the relic density, the coannihilation contribution has little to 
no effect on the relic density 85fl. In this region, the neutralino is dominantly Bino, making 
the coupling, and therefore, the coannihilation rate with the Wino-dominated X2 an d xf 
small 6 . 



sleptons. 

6 Again, assuming the sleptons are heavy, thereby suppressing the coannihilation process XiXi ~ * //■ 
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The calculated relic density of the neutralino in the various models is shown in Fig. [Tj 



versus the mass of Xi- The region 0.123 > QdmIt 2 > 0.099 is the 2cr observed range [72j, 173 ]. 
Values Q x oh 2 < 0.099 are allowed in the event that other particles contribute to the relic 
abundance of dark matter, or there is an enhanced non-thermal production mechanism. 

The profiles of the relic density versus the lightest neutralino mass are similar for the 
MSSM, NMSSM, and UMSSM due to their similar low energy neutralino spectra. In the 
NMSSM and UMSSM, the singlino is typically substantially heavier than the lightest neu- 
tralino, Xi! preventing a large singlino mixture in Xi- However, the NMSSM allows a lower 
neutralino mass than in the UMSSM and MSSM, as the k — ► limit resembles the n/sMSSM 
which has a very light Xi- 

Overall, there is a broad filled region in Fig. Q]for the MSSM, NMSSM and UMSSM where 
the relic density varies depending on the composition and mass of the lightest neutralino. For 
gaugino dark matter, the annihilation rate is not strong enough to yield the observed dark 
matter. As the higgsino content increases, the relic density decreases as the annihilation rate 
becomes larger. The suppression of the relic density due to efficient annihilation through 
the Hi pole near m x o ~ 60 GeV is evident. The H2 pole may also have a similar effect, as 
indicated by the sporadic points falling below the filled regions. In this low region, namely 
tt x oh 2 < few x 10" 3 , the NMSSM and UMSSM have more points than the MSSM since 
there are more Higgs resonances that contribute. 

Z X, X, Coupling 



mm?-. • 



30 40 
M (GeV) 



FIG. 2: The ZxiXi coupling (oc I-/V14I 2 — | A^i3 1 2 ) decreases in the n/sMSSM as the neutralino mass 
decreases, assisting in the enhancement of the relic density in the low mass region. 



For the n/sMSSM, the lightest neutralino mass is usually below 50 GeV 5l|, |80j. Anni 
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hilation through the Z boson pole dominates the rate near m x o « 45 GeV, decreasing the 
relic abundance below the observed value. Additionally, in the n/sMSSM the relic density is 
strongly dependent on the mass of the lightest neutralino. As the neutralino mass decreases, 
the annihilation rate becomes suppressed by the Z propagator, increasing the relic density. 
Furthermore, in the n/sMSSM, there is a Zx\x\ coupling suppression since the neutralino 
increasingly becomes singlino at lower m x o as shown in Fig. [2J Therefore, an approximate 
lower bound on the lightest neutralino mass in the n/sMSSM can be placed at m x o > 30 
GeV. 

Some points in the n/sMSSM exist below the 30 GeV m x o bound where the relic density 
falls at or below the observed value due to the contribution of additional channels. For 
example, a light A\ resonance enhances the annihilation cross-sections. Similarly, the A- 



funnel region in mSUGRA parameter space 86|, |87| can account for the observed relic density. 



However, in the present case, the lightest CP-odd Higgs is singlet. Excluding annihilation 
through the A-pole, the neutralino mass range compatible with the lower and upper bounds 
of Qdm are 30 GeV <; m^o < 37 GeV 7 . The limit might also be slightly weakened in the 



full secluded model 



35 



56| which may include a singlino Xo even lighter than the xi- Even 
a tiny Xo — Xi mixing, irrelevant for collider physics, could allow Xi to decay, e.g., to Xoffi 
reducing Q x oh 2 by m x o/m x o [881 ] . 

To further document the effect of the H\ pole on the relic density, we plot in Fig. [3] the 
masses of the lightest Higgs boson and neutralino for various ranges of the relic density. In 
models with singlet mixing, the lightest Higgs mass can be lower than in the MSSM, filling 
the major band of 114 GeV < M Hl < 135 GeV. 

In Fig. El the MSSM, NMSSM, and UMSSM show bands of allowed points associated 
with the Hi pole near the line m x o ~ 2.1mu 1 . The overdense region where m x o ~ 75 
GeV is due to inefficient annihilation of gaugino-like neutralinos. In the NMSSM, there are 
points with low allowed H\ masses near this band; most points there are due to annihilation 
through a SM-like H 2 . 

The n/sMSSM shows a weak dependence on the Higgs mass since the most dominant 
annihilation is through the Z boson. The vertical bands of accepted points are due to the 



7 Since we only consider sfermion masses much heavier than the neutralino, coannihilation with these 
particles do not contribute to the effective annihilation rate. In general slepton and squark coannihilation 
can be important as we parameterize the MSSM more generally than mSUGRA. 
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FIG. 3: Lightest Higgs mass versus lightest neutralino mass for the MSSM, NMSSM, n/sMSSM, 
and UMSSM. An over-density of relic neutralinos is denoted by a black x, an under-density is 
denoted by a blue box and the WMAP observed density is denoted by a red circle. 



dependence of the relic density on the lightest neutralino mass, as discussed above. 

Fig. H] shows the relic density in the plane of /i e fj and the lightest neutralino mass. Two 
bands of allowed relic density values are evident in the MSSM, NMSSM and UMSSM. The 
horizontal band at Mh ~ 60 GeV is associated with the lightest Higgs pole, while the other 



band is the focus point region 



89 



901 ] . This region is characterized by a Higgsino-bino mixed 



neutralino which makes neutralino annihilation efficient enough to reproduce the observed 
relic density |90J]. Note that even though the focus point region requires heavy sleptons 
and mixed higgsino dark matter, the neutralino may still coannihilate with neutralinos and 
charginos of similar mass. Increasing /i e g above this region makes the neutralino more 
gaugino-like, which decreases the annihilation rate. However, decreasing /j, e s will make relic 
neutralinos annihilate more efficiently, resulting in a relic density below the observed finish 2 
band. The n/sMSSM has a neutralino that is dominantly singlino, so the relic density is 
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FIG. 4: Regions of varying relic density in fj, e s vs. iti x o- The preferred horizontal strip in the 

» ~ 60 GeV is du 

region corresponds to the focus point. 



MSSM, NMSSM and UMSSM at m x o ~ 60 GeV is due to the Higgs pole. The steeply sloped 



insensitive to ji e s- 

The gaugino composition of the lightest neutralino is shown in Fig. where Z g = 
|iVn| 2 + | A^i2 1 2 is the MSSM gaugino fraction. The right sides of the panels indicate a large 
gaugino fraction while the left sides indicate neutralinos with high Higgsino or singlino/Z'- 
ino composition. In the MSSM, the broad band of WMAP allowed points between 1 < 
Z g / (1 — Zg) < 10 correspond to the focus point region while the other scattered points with 
higher gaugino fraction are due to the Higgs poles discussed earlier. 

In Fig. [6j we show the singlino and Z'-'m.o compositions of the neutralino. The UMSSM 
and n/sMSSM both have a range of singlmo/Z'-ino fractions that are dominant, but not 
close to maximal. As the neutralino becomes less MSSM-like, the relic density increases. 

Since the coupling between the neutralino pairs and SM particles (Z and Higgs bosons) 
needs a suitable value to give an acceptable annihilation rate, the lightest neutralino requires 
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FIG. 5: Neutralino relic density dependence on the MSSM gaugino fraction Z g . Points at the left 
end of each panel have a lightest neutralino that is dominantly higgsino and singlino, while at the 
right side it is dominantly gaugino. 

a nonzero MSSM fraction 8 . In the n/sMSSM, the annihilation is dominated by the Z boson, 
so the lightest neutralino requires a non- negligible Higgsino component. 

The singlino and Z'-ino dominated neutralino in the UMSSM can also explain for the 
observed relic density. The coupling between the lightest neutralino pair and the singlet 
dominated H 2 can be as strong as the XiXi^i coupling in the MSSM, c.f. Eq. (1B3[) . The 
strong coupling and resonant enhancement can yield a large enough annihilation rate to be 
below the Qdm observation. However, there are parameter points which do fit the observed 
relic abundance. 



Since we assume that the sleptons are heavy with masses at the TeV scale, the contributions of i-channel 
diagrams are suppressed. 
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FIG. 6: Neutralino relic density dependence on the singlino/Z'-ino fraction Z$- Points at the left 
end of each panel have a lightest neutralino that is MSSM-like, while at right side it is dominantly 
singlino or Z'-ino. 



III. 



XiP 



X?p ELASTIC SCATTERING 



The singlet extended models can give significant changes in the predicted MSSM cross- 
sections relevant to future recoil direct detection experiments. We use the modified version of 
DarkSUSY discussed in Section [III to calculate the spin- independent (scalar interaction) and 
spin-dependent (vector interaction) elastic scattering cross-sections of the lightest neutralino 
off nucleons. 

The experimental sensitivity to spin-independent (SI) scattering is much larger since 
spin-independent processes scatter coherently, and therefore are enhanced in scattering from 
large target nuclei. However, spin- dependent (SD) measurements can be made and have been 
probed with the Cryogenic Dark Matter Search (CDMS) to the few pb level for elastic proton 



scattering 



91] . The Chicagoland Observatory for Particle Physics (COUPP) experiment 
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which uses superheated CF 3 I can greatly improve the limits on the SD processes, down to 
10~ 2 pb for proton scattering for a 2kg chamber, which is close to the upper cross-section 



expected in the MSSM 

Current detection experiments such as EDELWEISS 



93] and CDMS have SI sensitivities 



on the order of 10 6 pb. In 2007 the sensitivity of the CDMS II experiment is expected to 



improve to nearly 10 8 pb 



941 ] . The proposed future upgrade, SuperCDMS, would reach a 



detection sensitivity of 10 9 pb. The WARP experiment with warm liquid Argon is projected 
;o reach a sensitivity of 10~ 10 pb and below; this experiment has just reported first results 



951 ] . The WARP sensitivity of an initial run with 2.3L and a total fudicial exposure of 96.5 



kg- day are slightly better than that o 



Dtained with EDELWEISS. For recent reviews on direct 



detection experiments see Refs. 



75 



m 



771 ] and for a forthcoming comprehensive summary 



of the status of direct detection experiments, see Ref. [191 ] . 



The SI scattering cross-section of a neutralino off a nucleus is given by 811 ] 



(i 



a% = ^\ZG° + (A - Z)G 



71 



n|2 



(4) 



where fj, x i 



m 0+m N 



is the reduced nucleon-neutralino mass. The parameters G p s and G" are 



hadronic matrix elements. Since the cross-sections for scattering off protons and neutrons 
are very similar, we focus on scattering from protons, for which the cross-section is 

,2 



a 



SI 

XP 



XP 
71 



\G p s \'' 



(5) 



Here G p s is given in terms of the hadronic matrix elements (N\qq\N) and couplings by 81] 

6 



G p s = £ (N\qq\N)-J2 

q=u,d,s,c,b,t k=l 



9q~LkXq9q~ Rk xq 



9 



(6) 



Ap) 

q h=H 1 ,H 2 ,H 3 



9hxx9hqq 



mi 



E 9h X x9hqq \ 2^ X - 
ml j 27 ^ 

h=H 1 ,H 2 ,H- i h / q=c,b,t 

where the three terms correspond to the diagrams in Fig. [7(a,b,c), respectively. We incor- 
oorated the following updated hadronic matrix elements to the DarkSUSY code taken from 



96] 



J: 



Tu 



(N\qq\N) 
0.020 ± 0.004, f v Td 



p n^p 



Tq 



m n 



TG 



0.026 ±0.005, f Ts = 0.118 ±0.062, 
0.84 



(7) 
(8) 
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(a) (b) (c) (d) (e) 

FIG. 7: Diagrams contributing to neutralino-hadron scattering. Scalar quark (a), Higgs boson (b) 
and Z boson (e) exchange contribute to the SI scattering cross-section. A significant contribution 
also comes from the Higgs scattering off gluons via heavy quark and squark loops (c,d). Processes 
which contribute to the SD scattering include squark exchange (a) and Z boson exchange (e). 
However, SD predictions are far less constrained than SI. 



We include effects from the exchanges of the scalar quarks of all six generations and the 
Higgs exchange from light quarks and from gluons via heavy quark loops. The diagrams are 
shown in Fig. [7J The scalar quark contributions are suppressed by our choice of TeV squark 
masses. The uncertainty in the SI scattering cross-section are large, of order 60%, due to 
the above uncertainties in the hadronic matrix elements. 

In Fig. [HI we show the predicted direct detection cross-section for the MSSM and the 
singlet extended models along with the sensitivities of EDELWEISS, CDMS II (2005), CDMS 
2007, SuperCDMS (25 kg) and WARP (2.3L). Over most of the region of neutralino mass, 
these experiments could find a signal for the MSSM, NMSSM and UMSSM. There are more 
points where the lightest neutralino of the MSSM may be detectable via recoil experiments 
than typical CMSSM expectations since the general parameter treatment of the MSSM can 



yie 



m 



d scattering cross-section that are a few orders of magnitude larger than in the CMSSM 



76 



rz|. 



Many of the points consistent with the observed relic density are within the reach of 
SuperCDMS and WARP (2.3L). The region below the WARP sensitivity corresponds to 
neutralino annihilation through the Hi (or H 2 ) pole. To account for the observed relic 
density, the XiXi^i coupling is small to balance the resonant enhancement of the annihilation 
rate. This forces the scattering cross-section via the Higgs states to be small, and thus fall 
below future sensitivities. 
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FIG. 8: Expected SI direct detection cross-section for (a) MSSM, (b) NMSSM (c) n/sMSSM 
and (d) UMSSM. The expected sensitivities of the EDELWEISS, CDMS II (2005), CDMS 2007, 
SuperCDMS (25 kg) and WARP (2.3L) experiments are shown. Over most of the neutralino mass 
range, the experiments should detect the signals from the MSSM, NMSSM and UMSSM. However, 
if the neutralino annihilates via a Higgs boson resonance, the relic density may be in the preferred 
region while the direct detection rate is out of reach of future experiments. 



Since the n/sMSSM has a light singlino state, the strict limits on the allowed neutralino 
mass in this model from the expected relic density can be used to deduce a lower limit 
on the direct detection cross-section. This is apparent in Fig. Ws) as most of the allowed 
region is above the expected sensitivity of CDMS 2007, which should allow CDMS to place 
extremely strong constraints on this model if no signal is found, assuming that the neutralino 
is responsible for a major part of the dark matter in the universe. Still, detection could be 
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precluded if annihilation through the Higgs resonance occurs with m x o < 30 GeV and 
Mffj ~ m x o/2, as shown by the isolated low Q x oh 2 points in Fig. [lb, or by neutralino decays 
into an additional still lighter particle [88]. 

The predictions of the SI cross-sections are heavily dependent on the local density of dark 
matter in our galaxy. The present halo model has many assumptions and uncertainties. 
Caustic rings are predicted to exist in the galaxy with a ring overlapping in the vicinity of 



our solar system 



97j . That could may increase the local dark matter velocity by a factor of 



three or more, resulting in a correspondingly enhanced flux and direct detection rate. 



IV. CONCLUSIONS 



In our investigation of the relic density constraints and the direct detection capabilities 
of present and upcoming dark matter experiments, we modified the DarkSUSY code to 
include the additional couplings and processes of sing let extended MSSM models 9 . We 
also updated the DarkSUSY values for the hadronic matrix elements and experimental 
constraints. Our analysis assumed multi-TeV masses of squarks and sleptons. We generally 
found close similarities of the NMSSM and UMSSM predictions with those of the MSSM due 
to the similarity of the compositions of their light neutralino states. However, the n/sMSSM 
showed exceptional differences due to very light neutralinos that are dominantly singlino. 
Some notable results from out study include: 

• The observed relic density can be accounted for in all of the singlet extended models, 
mainly with mixed Higgsino-gaugino dark matter. In some cases, the lightest neu- 
tralino is dominantly singlino, or a mixture of singlino and Z'-'mo in the case of the 
UMSSM. The predicted relic density can also match the observed value if the neu- 
tralinos annihilate through a Higgs boson resonance. The annihilation via the lightest 
Higgs boson can enhance the rate sufficiently to yield the correct relic density, even 
though the neutralino is typically bino in this case. 

• The neutralino in the n/sMSSM typically has a mass in the range 30 GeV < m x o < 
37 GeV to account for the total relic density found by the WMAP collaboration. The 

9 Lighter squarks and sleptons would allow the possibility of additional coannihilation regions. 
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upper bound on m x o may extend 



the observed DM value 



38, 



51 



;o 50 GeV if the neutralino relic density is below 
80] . Masses smaller than 30 GeV may be allowed if 
annihilation through the light CP-odd Higgs boson or neutralino decays occur to a 
still lighter state. 

• The singlino/Z'-ino dominated neutralino in the UMSSM can account for the observed 
relic density. Since the coupling between the lightest neutralino pair and the singlet 
dominated if 2 can be large and resonant enhancement of the annihilation cross-section 
via H2 can occur, the relic density can fall into the observed 2a range. 

• The MSSM, NMSSM, and UMSSM predict spin-independent proton scattering cross- 
sections that may be detectable at SuperCDMS and WARP and be consistent with 
the WMAP Qdm measurements. However, the recoil predictions of some models may 
be undetectable by these experiments due to the small Higgs neutralino coupling. 

• The n/sMSSM SI proton scattering cross-sections are highly favored to be detectable 
at CDMS 2007 while being compatible with WMAP Qdai observations. However, if 
neutralino annihilation occurs through a light Higgs, measurement of the scattering 
cross-section can fall below the sensitivities of future experiments. 

• Our MSSM predictions are more general than mSUGRA results. 

In addition to these constraints from the relic density bounds our study of further constraints 
in the Appendix found 

• Perturbativity constraints on A from RGE evolution give upper bounds on A(M t ) below 
unity and require tan/3 > 1.9. 

• Requiring that the xMSSM models are consistent with the anomalous magnetic mo- 
ment of the muon, Aa M , provides lower limits of tan/3 > 5, though there are still 
theoretical uncertainties. 
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APPENDIX A: NEUTRALINO MASS MATRIX 

The neutralino mass matrix of the singlet extended models is extended due to the addi- 
tional singlino and Z'-ino states that mix with the MSSM gauginos and higgsinos. In the 
(B°, W°, H Q d , H%, S°, Z'°) basis, the mass matrix is 

( Mt -giVx/2 giv 2 /2 \ 

M 2 g 2 v 1 /2 -g 2 v 2 /2 

-giv 2 /2 g 2 v 2 /2 -fi cS -fi c sv 2 /s g Z 'Q' Hl v i 

giv 2 /2 -g 2 v 2 /2 -// eff -n&V\/s g Z 'Q' H2 v 2 

~^eSV 2 /s -fieSVl/s \[2kS gz>Q's s 

o o gz'Q' Hl v i gz'Q' H2 v 2 gz>Q' s s M v 



V 



(Al) 



/ 



where the Mi are the gaugino masses (assumed universal at the GUT scale), and V\, v 2 and 
s are respectively y/2({H d ), (H u ), (S)). 

APPENDIX B: XMSSM COUPLINGS 

Couplings of the singlet extended models are changed with respect to the MSSM due to 
the additional singlet and Z' contributions. 



The Yukawa the CP-even and CP-odd Higgs couplings are 



- n SM 

gddH, — 9/fh 



R 



il 



gddA, = «7 9fth- 



cos (5 
R 1 } 



- n SM 

guuHi — g/fh 



R 



■42 



sin (3 



guuA, = rfg s f fl- 

J J Q-l 



R 



i'2 



:bi) 



(B2) 



cos/? ^ ' ^ ft sin/?' 
where Rf. and R % i are the rotation matries that diagonalize the CP-even and CP-odd 
Higgs mass-squared matrices, respectively, in the (H u , H d , S) basis, and g§M are the 
corresponding couplings in the SM (see Ref. [371] ). 



There are some couplings that do not have a standard model counterpart, as follows: 
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The Higgs-Neutralino-Neutralino coupling constants are: 

9H tX Wi = [(9iNii-g2N 12 -g v Q Hd N m )N 13 + V2XN 14 N 15 \R i l 

+ [Mfo - gxNxx - gvQ Hu N 16 )N u + V2XN 13 N 15 ]R^ (B3) 
+ [-gvQ S N 16 N 15 + V2\N 13 N U - V2kN 15 N 15 ]R™ 

9A iX U°i = rf[(9iN 1 i-g2N 12 -g 1 ,Q Hd N 16 )N 13 + V2\N u N 15 ]R 1 } 

+ [(g2N 12 - 9l N u - grQ Hu N 16 )N 14 + V2\N 13 N lb ]R^ (B4) 
+ [-9vQ S Ni 6 N ls + V2XN 13 N 14 - V2kN 15 N 15 }R*, 

where Nij is the rotation matric that diagonalizes the neutralino mass matrix in Eq. 

The Left and Right Handed Neutralino-Quark-Squark coupling constants are: 

gLu kX o Ui = —^(-N n9l + N 12 g 2 ) - l^zf 3)fe *iV 14 (B5) 

gm kX « Ul = ^zt 3)k *N* u - Y^tm, (B6) 
Z ik 1 

9Ld M = -^(-N ll9l - N 12 g 2 ) - Y*zl +3)k N 13 (B7) 



gm M = -^Z^ k N* n - YiZfN{ 3 (B8) 



where Z~ and Z~ are the rotation matrices for the left and right handed scalar 
quarks, respectively. The Yukawa couplings, Y*, are defined for quark qi. 

The trilinear Higgs coupling is modified due to the additional singlet state in the 
nMSSA/LNMSSM, and UMSSM. The HiAjA k coupling can be found using the pro- 



jection 



37] 



C HiAj A k = P Hi P Aj P Ak V, (B9) 
where V is the Higgs potential and the projection operators are 

p ^U R "^ +Rt2 ik +Ri 4) (Bn) 
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In a similar manner the HiHjH^ coupling is 

C HiHjHk = P Hl P Hj PH k V, (B12) 



and the HiHjAk coupling is 



C HiHjAk = P Hi P Hj PA k V. (B13) 



APPENDIX C: EXPERIMENTAL CONSTRAINTS 

We generated the Higgs and neutralino masses and couplings by scanning over the pa- 
rameters given in Ref. [371] and diagonalizing the Higgs and neutralino mass matrices. For 
each model, the following LEP experimental constraints are applied: The ZZh coupling 
limits, the bound on the lightest chargino, the charged Higgs mass bound, the Ah associated 
production search limits and the contribution to the invisible decay width of the Z boson, 
and the limit on the Z — Z' mixing in the UMSSM, as in Ref. [371] . In addition, we con- 
sidered limits from the anomalous magnetic moment of the muon, coupling perturbativity 
constraints and the improved limits on BF(6 — > 57) [98 ]. 



1. Anomalous Magnetic Moment of the Muon 



The experimentally measured value o: 
= (9 



2)^ shows a 3.4cr deviation 



99 



the anomalous magnetic moment of the muon 



100| 



Aa M = cv(exp) - a At (SM) = (27.6 ± 8.1) x 10 



-10 



(CI) 



from the Standard Model prediction based on the use of e + e~ — > tttt data from CMD-2 
and KLOE. The constraint s on the UMSSM parameter space due to the a M deviation were 
previously studied in Ref. 1011 ] . We apply the above a M constraint here and find that the 
results do not vary significantly across models. In particular, Aa^ gives modest lower limits 
on tan (3 when using the parameter values in Ref. [37J in a scan over the soft slepton masses 
300 GeV < M~ L = M E < 2 TeV. At the 2a level the Aa M bound limits tan/? > 5. We did 
not impose this constraint in our analyses, because there are still theoretical u ncert ainties 
involving the disagreement between estimates using e + e~ and r decay data 
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FIG. 9: The deviation of the anomalous magnetic moment of the muon Aa^ consistent with the 
2<T uncertainty for values of tan (3 in various models is within the dashed horizontal lines. 

2. Perturbativity constraints 



The singlet models that we consider may have gauge coupling unification at a Grand 
Unification scale Mqvt ~ 2 x 10 16 GeV. We may require that the gauge and Yukawa 
couplings remain perturbative up to this scale. To impose this constraint, we evaluate the 
running of the cou plings in each model according to the re normalization group equations 
(RGEs) 



46 



58 



59 



60 



74 



102 



103, 



104 



105 



106 



107 



10 a ] at one loop order: 



dg% 
dt 
dgv 
dt 
dn 

~dt 
dX 

dh t 
~dt 



Pi9i, 



6k [X 2 + k 2 ] 



A 

h t 



4A 2 + 2k 2 + 3h 2 
16 



3gl 



^gf-ZgUQl + Qk + Qk 

13 



A + Q 1 ~ Y 93 ~ 392 ~ TE 91 ~ 2gi ' iQu + Qq + Qh - 



(C2) 
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where t = j^K? ln(fi/m t ). Here h t is the top quark Yukawa coupling and Qh u , Qii d , Qs, Qui 
and Qq are the U(l)' charges of the up and down Higgs doublets, the Higgs singlet, the up 
type quark and the quark doublet, respectively. For specific models, the parameters g[ and 
K are appropriately turned off. 

The scale factors and /3-f unctions are given in Table II VI 10 . In Table II I II we give the 
quantum numbers under the Stand ard M odel gauge groups, U{l) x and U{1)^ (symmetries 



arising when the Eq model is broken 1021 ]). and of U(l)' (the additional symmetry of interest 
in the UMSSM). As U(l)' is simply a linear combination of U(l) x and U(l)^, the quantum 
numbers under this symmetry are given by 

Q = Q x cos6 E 6 + Qipsin6 Ee (C3) 

and thus depend on the particular value of 9e6 ■ 

The RGEs were evolved from the low scale, chosen to be m t = 171 GeV, up to Mqut ~ 
2 x 10 16 GeV. The gauge couplings were first run independently up to m t from Mz, using 



22, 



HQ I 



sm 2 9 w (M z ) = 0.23122 ±0.00015 

a(M z ) = 1/(127.918 ±0.018) (C4) 
a s (M z ) = 0.1189 ±0.0010 

where a = e 2 /4n, as = c?|/47r, g x = \j\gy = \f\ e l cos6 W} and g 2 = ej sin^^. The 
gauge couplings were assumed to unify at the GUT scale, with the canonical normalization 
9v — 9i — 92 = g-i at the unification scale 11 . From the top quark pole mass, we obtained 
the low-scale value of the top quark Yukawa coupling 

7 / x rn t (m t ) / r ,, 
h t {m t ) = — C5 

For the couplings to remain perturbative up to M GUT , we require 

k 2 < 4vr ; A 2 < 4vr ; h 2 < 4vr (C6) 



10 We do not include the effect of kinetic mixing in the UMSSM 

11 0y d epends weakly on 8e 6 , which is scanned, because of the H 3 ,H 3 pair needed for gauge unification 



102[. This implies that gy = •jXggx at the weak scale, where X g , which depends on 9e 6 , is typically 



within 6% of unity. 
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FIG. 10: Values of A vs. tan/3 consistent with all considered constraints except perturbativity in 
scans over all parameters. The solid and dashed lines represent the maximum allowed values of 
X(rrit) that remains perturbative in the NMSSM for k = and k = 0.5. The n/sMSSM constraint 
is given by the k = result of the NMSSM. 

The NMSSM perturbativity constraints are shown in Figure \IU\ as the solid line for k = 
and a dashed line for n = 0.5. The n/sMSSM corresponds to the k = limit of the NMSSM. 
There are many points with large values of A at low tan f3 that violate the limit in the k = 
NMSSM and n/sMSSM. Figure [TUl also demonstrates that perturbativity becomes slightly 
more constraining in the NMSSM for larg QT Kj^ clS Ki contributes to the running of A. The 
maximum allowed values of X(m t ) for moderate values of tan/3 are given in Table HT1 for 
several values of k. These values are slightly lower than the constraints of A(m t ) < 0.75 and 
^\ 2 {m t ) + k 2 < 0.75 for the NMSSM that are often used as approximations to the RGE 
running |37l] 12 . 

The perturbativity limit for the UMSSM is given in Figure [IT] for the example 13 of 



12 Note that these constraints can be weakened if one adds ad ditional vector- like matter, increasing the 
(negative) effect of the gauge couplings in IC3I 
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The value for which the U(l)' charge of the right handed neutrino Nf vanishes. 
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FIG. 11: Values of A vs. tan/3 consistent with all considered constraints except perturbativity in 
scans over all parameters in the UMSSM, including 9 eg- The solid line represents the maximum 
allowed values of A(mt) that remains perturbative for the particular choice of 9 eg = arctan\/l5. 



321, 



Unlike the other models, A is strongly constrained by the requirement of a small Z — Z' mixing 
which, for reasonable /^ e fj, favors large {S) and small A. This is especially true for large tan/? or 
small 9e 6 , for which there are no cancellations between the contributions of H% and H® to the 
mixing. 



TABLE II: The maximum allowable value of A for perturbativity in the NMSSM, as a function of 

K. 



K 


0.0 


0.3 


0.4 


0.5 


0.6 


Amax 


0.70 


0.64 


0.58 


0.47 


0.25 



9eg — arctanvl5. It is less constraining than the NMSSM, with very few models being 
eliminated due to A moa: (mt), which varies from ~ 0.82 to ~ 0.86, depending on the choice of 

In addition to the perturbativity upper bound on X(m t ), the RGE running sets a lower 
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TAB LE III: Matter multiplets in the UMSSM and their various quantum numbers (according to 
1021 ]. but noting our difference in normalization). The index i sums over generations. 



super fields 


SU(3) C 


SU(2) L 


Qy = \f\Q\ 


17(1)' 


\/40Q x 


V^Q^ 


4 


3 


2 


1 

6 


Qq 


-1 


l 


n 


3 


1 


2 
3 


Qu 


-1 


l 


k 


3 


1 


1 

3 


Q d 


3 


l 


u 


1 


2 


1 
2 


Ql 


3 


l 


k 


1 


1 


1 


Qe 


-1 


l 




1 


2 


1 
2 


Q Hi 


-2 


-2 




1 


2 


1 

2 


Qh u 


2 


-2 


Di 


3 


1 


1 

3 


Qd 


2 


-2 


A 


3 


1 


1 

3 


Qd 


-2 


-2 


Si 


1 


1 





Qs 





4 


Nf 


1 


1 





Qn 


-5 


1 


H 3 


1 


2 


1 

2 


-Qh u 


-2 


2 


H 3 


1 


2 


1 

2 


Qh u 


2 


-2 



TABLE IV: One-loop (3 functions for the three models 102l . 1 1031 ] . The traces involve sums over all 
matter superfields and generations. 





0i 






$v 


NMSSM/nMSSM 


33/5 


1 


-3 





UMSSM 


iTrQl = 9 + f 


4 





TrQ 2 = 9 + 4Q^ 



limit of tan > 1.8-1.9 for the NMSSM/nMSSM and tan/5 > 0.8-0.9 for the UMSSM. Be- 



low these tan (3 values, a large h t often feeds into the runn ing o 



and sometimes resulting in a Landau pole below M GUT 



103 



A, making A n on-p er t urbat ive 



105, 



106 



107]. 
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